Thermal hadron spectrum in e + e annihilation from gauge/string duality 
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behavior. 
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Introduction — One of the puzzling features of hadron 
production in collider experiments is that the inclusive 
rates are well described by assuming the 'Boltzmann' dis- 
tribution at low momenta 
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where the subscript i labels different hadron species. 
The parameter T, often referred to as the 'tempera- 
ture', is more or less independent of the collision en- 
ergy and is typically of the order of the confinement scale 
> 150 MeV. The model (JTJ gives a good fit of the iden- 
tified particle yields in e + e~ annihilation [l|, 0, Q (see, 
however, Q), hadron collision [j| and heavy-ion collision 
experiments [1, 0] ■ Despite these phenomenological suc- 
cesses, however, the origin of the exponential behavior is 
not understood. One may invoke Hagedorn's picture of 
particle production Q in which an equilibrated state of 
secondaries ('fireball') is instantly formed after the col- 
lision. Although such a scenario may sound plausible 
in the heavy-ion case where one expects the formation 
of the quark-gluon plasma, the situation is highly ob- 
scure in the case of e + e _ or hadron collisions where the 
produced particles have essentially no chance to interact 
among themselves. 

In this letter we derive the exponential distribution 
in e + e~ annihilation from gauge/string duality, thereby 
suggesting that the apparent thermal behavior is a 
generic feature of the strong coupling dynamics of gauge 
theories in the nonperturbative regime. [See [9( for 
a quite different approach to par ticle production using 
gauge/string duality. See also jiol . Ull . Il2l . [l3| for earlier 
speculations concerning the behavior (|T]).] e + e~ anni- 
hilation in J\f = 4 supersymmetric Yang-Mills (SYM) 
theory and its variants has been recently analyzed in 
the framework of the AdS/CFT correspondence [HI]. 
Here we briefly summarize the key features at large 't 
Hooft coupling A 3> 1. (i) Unlike in QCD and in 
J\f = 4 SYM at weak coupling, there are no jets in 
the final state [HI, [HI, [T3|- The distribution of energy 
is spherical [HI], (ii) The average multiplicity scales al- 
most linearly with the virtual photon energy Q 0, [TlJ : 
n(Q) tx (Q/A) 1 " 3 / 2 ^ where A is an infrared cutoff, (iii) 
The inclusive spectrum is peaked at the kinematical lower 



limit. More precisely, it has the form [l7| 
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where x = 2E/Q is the Feynman variable, 
distribution is spherical, this is equivalent to 
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Since the 
(3) 



The function F(y) remained unknown except for the 
property that it decays faster than the power-law when 
y 3> 1. In 17], it was conjectured that 
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where c is a number of order unity. We shall show that 
a calculation based on gauge/string duality yields the 
value a = 1 in accordance with the phenomenological 
distribution (T]) with 1/T = 2c/A. 

Statistical model — As a preliminary, let us first point 
out the possible relationship between particle production 
in strongly coupled gauge theories and a statistical model 
of e + e~ annihilation proposed by Bjorken and Brodsky 
long ago [lj[. Suppose, as is the case in QCD, that the 
multiplicity is dominated by the lightest particles ('pi- 
ons') with mass m. The inclusive distribution is given 
by 
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where k and k! and electron and positron momenta, re- 
spectively and jfj, is a component of the 7?.-current opera- 
tor which couples to the external virtual photon. Rather 
than directly dealing with this expression, following [l8| 
we consider the cross section of producing exactly n pions 
and relate it to the inclusive distribution. 
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We take n to be large, of the order of the average mul- 
tiplicity n ~ Q/A » 1. The model assumes that the 
hadronic matrix element can be written as 



(0|j"(0)|pi,. ..p„)(pu- Pn\f (0)\0) 
-a n ( g V-5' 1 V)e^ (Bl+ - +E " ) , 
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where j3 — 1/T and the arrow means that the replace- 
ment is valid only under the n-particle phase space inte- 
gral. Note that the matrix element itself gives the 'ther- 
mal' factor, which distinguishes this model from other 
statistical models where the factor comes from the total 
phase space. [See also, pjj].] We then get 
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Due to the non-renormalization theorem for the 1Z- 
current correlator [20| , the total cross section is given 
by the one-loop result for all values of the coupling 

e 4 7V 2 

a tot — nn ■ (8) 
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After an integration over Q within an interval AQ <C Q, 
([9]) takes the form 
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where we approximated E' «p'. The last factor may be 
omitted and the summation simply gives a factor An ~ 
AQ/A 
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The coefficients {a„} can be determined by matching 
with the total cross section ([7]) [l8[ 
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[We have used Stirling's formula nl « v / 27rn n+1 / 2 e n 
valid for n>l.[ Inserting this into (fll)]) . we obtain 



d 3 p 7rA 7rA 2 
This agrees with ([3]) provided a — 1 and 
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The proportionality constant is determined from the en- 
ergy conservation 
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The result is 
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and for the average multiplicity, 
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It is straightforward to include heavier hadrons with mass 
to* > to as a small contamination in the final state. Pro- 
ceeding as before, we find 
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[Since typically I <C n, we assume that the /-dependence 
of a n is negligible, or trivially factorized in the form 
a n,i ~ with 6 being a number of order unity] This 
immediately leads to the characteristic ratio 



( n *) = m *K\ {(3m*) ^ e _ /3(m ._ m) 
(n) mi^i (/3to) 
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which resembles the observed scaling [T|, H, H, H, 0]- 
Therefore, finding the exponential factor in (J5]) 

establishes a connection between strongly coupled N = 4 
SYM and a statistical model of particle production in 
e + e~ annihilation. 
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Gauge/string duality — Let us now turn to the matrix 
element 



(0|e-i(0)|pi,...p n ), 



(19) 



where is an arbitrary polarization vector which we 
take to be transverse (e ■ q — 0) for convenience. The 
prescription for evaluating amplitudes like (fT§)) has been 
spelled out in [2~ll. [2^|. In the limit A — * oo, this can be 
viewed as the scattering of n + 1 bulk fields in AdSs x 
S 5 in the supergravity approximation. We shall use the 
Poincare coordinates 



ds 2 = R 2 
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The Minkowski boundary is at z — and the space is 
cut off at z = 1/A. The current j M is dual to a five- 
dimensional Kaluza-Klein photon whose wavefunction 
solves the Maxwell equation in the bulk. For timelike 
momenta, q 2 = Q 2 > 0, and with the infalling-wave con- 
dition as appropriate for the problem of jet evolution 
it features the Hankel function 
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[The normalization is taken from [23[ . Note that A z = 
in the current gauge (e • q = 0).] Alternatively (j2T|) 
may also be regarded as the wavefunction of a heavy 
vector meson with which the photon couples in a way 
similar to the vector meson dominance in QCD. Then the 
factor N c in ([2Tjl acquires the meaning as the strength 
of this coupling. Naturally, the timelike photon or the 
meson decays, and at strong coupling the decay is so 
complete that the final state contains only particles with 
the smallest virtuality p 2 = m 2 ~ A 2 [H, 0| . These are 
represented, for simplicity, by the scalar wavefunction 
in AdSs. We shall consider two scenarios, (i) Hadrons 
corresponding to normalizable modes in the cutoff AdS 
space. They have the wavefunction 
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where A + is the conformal dimension of the interpolating 
operator. The lightest one is dual to the lowest dimension 
scalar operator in the TV = 4 graviton supcrmultiplet. 
This has A + = 2 and m w 2.4A as determined from 
the first zero of the Bessel function Jo. (ii) The non- 
normalizable mode with A_ = 4 — A + = 1 
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This somewhat peculiar choice is motivated by the follow- 
ing argument. As noted in [24[ and emphasized in [25| . 
the value A_ = 1 realizes Bjorkcn scaling in DIS which 
is the hallmark of the parton picture. [See, also, [26j].] 
Curiously, these modes contain the so-called singletons, 
or rather, doubletons [23 that are often referred to as the 



most elementary excitations, or 'partons' in AdS spaces. 
Moreover, the bulk supergravity modes such as ([2"Tj) are 
bound states of singletons [H, [2!| . This fits nicely with 
the branching picture at strong coupling proposed in [l6| . 
Note that at the borderline value A_ = 1 the normal- 
ization integral suffers a logarithmic UV (z — > 0) diver- 
gence. This is analogous to the fact that the light-cone 
wavefunction with soft partons (gluons) in QCD is not 
normalizable due to the same reason [30j . 

In the branching picture of |l6j], one evaluates diagrams 
with only the three-point (1 — > 2) vertices. On the other 
hand, the supergravity effective action contains all possi- 
ble higher-point vertices allowed by symmetry. To better 
illustrate the essential point of our calculation, let us first 
consider a contribution to the matrix element from a lo- 
cal (n + l)-point vertex of n scalars and a photon. The 
zero mode integration reads, symbolically, 
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where g c ~ g s ct' 2 is the closed string coupling. The sim- 
plest vertex is obtained by expanding the function F to 
linear order in its argument 
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where v a (a — 1, 2, .., 5) is a Killing vector on S 5 and Qij 
are the U(l) charges. The sum is over oppositely charged 
pairs. Since n ~ Q/A is large, one can do the saddle 
point approximation. It is consistent to look for a saddle 
point in the asymptotic region of the Hankel function 
H^\Qz) ~ e l Q z / ^Qz. The equation determining the 
saddle point is, setting n = kQ/m where A: is a number 
of order unity, 



— — In (z 2 J (mz)) 
m dz 
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and similarly for H\ (Az) in which case we set n — 
kQ/A. We solve this equation numerically in the form 

a 

z=jt, 

and find that, in the physically interesting range 0.5 < 
k < 1.0, a w 0.2 - (0.8 ± 0.1)i and a « 0.5 ± 0.1 for 
the Jo and cases, respectively. This leads to, using 
Ang s a' 2 /R A = 1/N C and Oi = -Q h 
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where the factor 1 / y/Q/z is from the fluctuation around 
the saddle point. Squaring and noting that p%Pj — > 
g^Sijpf ~ g^SijA 2 under the phase space integral, we 
find 



\(0\e-j\ Pl ,...p n )f 
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where c = Re [a] . (|2T[) features the exponential factor as 
in ([5]). Note that having n ~ Q particles in the final state 
is crucial in order to obtain this behavior. If n is order 
unity, the integral would be dominated by intermediate 
z values, leading only to a power-law (2l| . 

The evaluation of the nonlocal diagrams with the 1 — > 
2 vertices is considerably more complicated, but without 
going into the details one can argue that they also give 
rise to the exponential factor. [The exact value of c may 
be different, however.] Indeed, as one approaches the root 
of the branching tree, the intermediate states necessarily 
have an increasingly large number of constituents, hence 
are dual to operators with increasingly large dimensions 
whose wave functions are localized around z ~ 1/A. In 
particular, the photon couples to fields &x.y having di- 
mensions of order Ax,y ~ n ~ Q/A via the vertex 
A^i^xd^Y-^vd^x) [See [H[ for a related discussion 
including the issue of gauge invariance of this coupling.] 
The matrix element will be given by a multiple convolu- 
tion of the form 

{Px-Py) j dzdz'dz" ■■■V^GA i1 (z)Gx(z,z')Gy(z,z") 

■ ■ ■ G(z"', ziMziMzi) ■ ■ ■ G(z"", z n/2 )$(z n/2 )$(z n/2 ) , 

where G's are propagators in (cutoff) AdSs. At small z, 
Gx,y(z, z 1 ) ~ z Ax - Y <~ z n , so the function obtained after 
the integration dz'dz" ■ ■ ■ is as strongly peaked around 
z~ 1/A as $ n (z) in (|24|) . The final integral over z then 
leads to the exponential factor by the same mechanism 
as in the local case above. 

In fact, although the two contributions-local and non- 
local diagrams-have similar Q-dependence, they cru- 
cially differ in iV c -dependence. When comparing © and 
(fTTjl with (|27p . one notices that the latter is enormously 
suppressed by the factor l/N 2n ~ 2 (1/Af" _1 in the am- 
plitude) , which actually was expected from the standard 



l&rge-N c counting. If the non-local diagrams also re- 
ceived this suppression, one would have to conclude that 
the multi-particle production process would be a negli- 
gible fraction of the total cross section |(7J) which scales 
as N 2 . Instead, the production cross section of a single 
heavy vector meson, being 0{N 2 ), would alone seem to 
saturate the total cross section. In order to avoid this 
unphysical conclusion, the subsequent decay of the me- 
son into n particles via 1 — > 2 splittings must somehow 
be an 0(N°) effect rather than suppressed by a large 
negative power of N c . Fortunately, a similar problem 
was encountered and resolved in the context of e + e _ an- 
nihilation in (two-dimensional) QCD at large N c [3l| . 
There it was observed that a proper inclusion of the de- 
cay width (r ~ 1/N 2 in the present situation) in the 
intermediate propagators together with a suitable aver- 
aging of the incident energy give just the right number 
of N c factors to cancel the unwanted 1/N C suppressions 
from the three-point couplings g ~ l/N c . Namely, the 
squared propagator in the momentum space reads 
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(p 2 — to 2 ) 2 + r 2 r 
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The first factor on the right hand side is unit normalized 
with respect to the integration over the momentum and 
the second factor is 0(N°), so effectively (2g) is 0(N°) 
contrary to the naive expectation 0(1/N 2 ). Since this 
argument appears to be quite general, we expect that 
the nonlocal diagrams will not be A^-suppresscd unlike 
the local ones. 

Keeping this qualification in mind, we conclude that 
the distribution obeys the thermal law |T]), and espe- 
cially, ([IH]), though clearly the parameter T ~ A has 
nothing to do with the temperature. In view of the strik- 
ingly different pictures of the final state in QCD and in 
strongly coupled SYM as mentioned in the introduction, 
it is intriguing that the latter can explain this particular 
nonperturbative aspect of particle production. 
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